Phonon Softening and Dispersion in EuTiOs 
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We measured phonon dispersion in single crystal EuTiOs using inelastic x-ray scattering. A 
structural transition to an antiferrodistortive phase was found at a critical temperature To=287=bl 
K using powder and single-crystal x-ray diffraction. Clear softening of the zone boundary i?-point 
q=(0.5 0.5 0.5) acoustic phonon shows this to be a displacive transition. The mode energy plotted 
against reduced temperature could be seen to nearly overlap that of SrTiOs, suggesting a universal 
scaling relation. Phonon dispersion was measured along F-X (0 0)^(0.5 0). Mode eigenvectors 
were obtained from a shell model consistent with the q-dependence of intensity and energy, which 
also showed that the dispersion is nominally the same as in SrTiOs at room temperature, but 
corrected for mass. The lowest energy optical mode, determined to be of Slater character, softens 
approximately linearly with temperature until the 70-100 K range where the softening stops, and 
at low temperature, the mode disperses linearly near the zone center. 

PACS numbers: 63.20.dd, 64.60.Ht, 77.80.e 



Materials in the ^TiOa family of perovskites, where 
A = Sr, Ba, Pb and others, have been studied with in- 
elastic neutron scattering [1-4] , with particular attention 
to soft phonon modes associated with structural and fer- 
roelectric transitions [sl, 0] . Although known for nearly 
60 years [7], EuTiOa had not been included in similar 
studies, in part due to the neutron-absorbing Eu atoms 
making such experiments difficult, and to the unavail- 
ability of bulk single domain samples necessary to resolve 
phonons in momentum space. It is becoming increasingly 
apparent, however, that EuTiOa has fascinating proper- 
ties, owing to the large magnetic spin moments on the 
Eu^+ ion which can couple to the electrical polarization 
This makes knowing the lattice dynamical prop- 
erties of this material all the more imperative, but studies 
have been largely limited to theoretical efforts [l3l - [l6| 
and, on the experimental side, zone-center infra-red re- 
flectivity [H, 1S\ which probes zero momentum only. The 
general consensus is a tendency towards ferroelectricity 
governed by softening of a zone center F point optical 
phonon, which was verified with the infrared reflectivity 
measurement. 

Recently, a peak was detected by Bussmann-Holder et 
at. in the specific heat of a ceramic sample of EuTiOa, 
at r=282 K [14]. They hypothesized it originates from 
an antiferrodistortive transition to a tetragonal 14/ mem 
phase. This transition features a quasi-rotation of the 
oxygen octahedra about one of the axes and often a c/a 
lattice constant distortion, and in SrTiOa exhibits an as- 
sociated soft phonon at the zone boundary R point in 
momentum space [l^, Since there have been only 

a handful of materials which clearly showed softening 
phonons related to this antiferrodistortive transition, and 
underdamped cases besides SrTiOs have been quite rare. 



EuTiOs therefore potentially represents an important ad- 
ditional window on the dynamics of such soft mode phase 
transitions. However, the analysis of [14], which correctly 
predicted the higher transition temperature for EuTiOs 
due to spin-phonon coupling, at the same time suggest 
that the transitions should have an order-disorder char- 
acter, as compared to the displacive transition of SrTiOs- 
Subsequent x-ray diffraction measurements of EuTiOs 
powder samples conffimed the 14/ mem low temperature 
phase for EuTiOs, but only below 235 K The au- 

thors of that study argued from pair distribution function 
analysis that the distortion initially forms in nm-size do- 
mains which only order at the lower temperature. How- 
ever, a subsequent study [27] observed a transition with 
powder measurements at around r=300K. 

To study its dynamics, we grew a bulk single crystal of 
EuTiOs for inelastic x-ray scattering, a technique which 
both bypasses the problem of neutron absorption, and 
relaxes the need for a large sample. Our single crystal 
samples were grown by the traveling float-zone method, 
using a variant of the method of Katsufuji and Tokura 
[221] . Full details of the growth and characterization 
will be published elsewhere, but a summary is given 
here. Feedrods were pressed from a mixture of the 
component powders Eu20s and Ti02, and sintered to 
form Eu2Ti207. Growth was in a reducing 95:5 Ar:H2 
gas mixture, using a commercially obtained SrTiOs 
seed crystal for the initial run, and iteratively using the 
top part of previous runs as seeds for the next run. A 
SQUID magnometer was used to measure the magnetic 
susceptibility of the resultant crystal, which showed a 
peak at Tat =5. 7 K, in accordance with the known Neel 
temperature of the antiferromagnetic transition [23]. A 
cleaved surface resulted in the best mosaic, typically 
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FIG. 1: (Color Online) (a)-(f) Selected phonon spectra at 
the i?-point (Q=(4.5 0.5 1.5), with resolution ^Q=(0.03 0.02 
0.03)) as a function of temperature as indicated. The solid 
lines are fits to a phonon mode (with no central component) 
as described in the text. 



0.01-0.03°. We used a (1 0) crystal face, which was 
verified to be consistently single domain over the surface 
of ~0.5 mm dimension. At temperatures below r=250 
K, the crystal became obviously twinned, but the overall 
mosaic (< 0.2°) and splitting in 20 (|AQ| -0.03 A ) was 
deemed sufficiently small to neglect for this study. 

After observing the twinning, we performed powder 
diffraction at BL44B2 of SPring-8 [24]. The results 
were broadly consistent with those in [21| . Our Rietveld 
refinement [25|, showed a transformation from a cubic 
phase at room temperature to a tetragonal 14/ mem 
phase below 200 K [26|. Broadening of expected twin 
peaks was visible below —240 K as in [21], however 
weak i^-point (1.5 0.5 0.5) intensity was visible in 
our powder data below —280 K, whereas it was only 
seen at —200 K in [21]. We also investigated the 
i?-point (3.5 0.5 0.5) reflection using single crystal 
diffraction. The intensity had a clear onset at 287±1 
K (see inset to fig. [2j), which is within a few degrees 
of the previous calorimetric data [1^. The angular 
width of the diffracted beam was largely temperature 
independent, varying from 0.23 to 0.33 mrad depending 
on location on the sample, at both 280 K and 220 K. 
From this we estimate a minimum correlation length 
>0.1 jiia. This suggests that the domain sizes are 
large, which does not appear consistent with disorder 
of the oxygen octahedra rotation angles suggested in |21|] . 
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Inelastic x-ray scattering measurements were done 
using BL35XU at SPring-8 |28|, using a 21.7 keV 
incident energy with typically 1.5 meV resolution and 



FIG. 2: (Color Online) Energy of the soft mode at Q=(4.5 
0.5 1.5), fitting the spectra using two methods. The empty 
squares are from fitting a single mode (a stokes/anti-stokes 
pair), and the filled circles correspond to including a central 
resolution- limited peak. As shown in the inset, To was deter- 
mined from the Bragg peak intensity as a function of tem- 
perature measured at the Q=(3.5 0.5 0.5) i^-point. The con- 
tinuous curve corresponds to SrTiOa [19'], for which To=108 
K. 



momentum resolution of 0.03, in units of reciprocal 
lattice. The i?-point spectrum at Q=(4.5 0.5 1.5) was 
measured between r=370 K and T=289 K at intervals 
of 10 K. Selected spectra are shown in Fig. [T]^a)-(f), 
which clearly show a softening mode. The data were 
fitted to two Lorentzians, whose peak intensities are 
constrained by detailed balance, and convolved with 
the resolution function. The fitted energies are plotted 
against reduced temperature t= ^~^^ in Fig. [21 In fitting 
the spectra we first neglected, and then considered, a 
possible "central mode" at zero energy [29|, (SOj. The 
mode energies resulting from fits which include it as 
a resolution-limited peak of freely varying intensity, 
are also plotted in Fig. O Near the transition, the 
energy derived without including a central peak is linear 
with temperature, and fitting to a line gives To=287 
K. This is consistent with value of To=287ztl K as 
determined from the Bragg intensity shown in the inset 
of fig. [21 For a comparison with SrTiOs, we plot its 
soft mode energy in fig. [21 using the empirical formula 
of Shirane and Yamada ^], who determined To=108 
K. With the second fitting method, the approach to 
To in the two materials becomes identical within error 
bar for ^^'^^ <0.1. Previously, Miiller and Ber linger 
had found a universality in the critical exponent of 
the order parameter below the transition temperature, 
by comparing SrTiOa to LaAlOa [31], but the present 
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result, assuming the plausible existence of a central 
component, suggests a further universality in the soft 
mode energy vs. reduced temperature. That the 
absolute phonon energies would line up may stem from 
a similarity between the two systems. 
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FIG. 3: (Color Online) (a) (bottom panel) Phonon energy vs. 
momentum q along the T-X direction. The dispersion away 
from the T-point (g>0.1) was for (H K L) = {4, 0, 1-g). For 
the near-F dispersion, zones (4 1-1) (empty triangles and di- 
amonds) and (5-1 -2) (filled triangles) were used as indicated 
in the legend. The smooth curves are model calculations, with 
the thickness of the curve proportional to the square root of 
the intensity in the (4 1) zone. The inset shows the spectra 
at the anti-crossing near 20 meV and g=0.25 at T=10 K. For 
clarity, the two distinct mode energies at g=0.26 have been 
slightly offset along the q axis, (b) (top panel) Relative inten- 
sities of the transverse acoustic and high-energy modes along 
(4, 0, 1-g), comparing the experimental values obtained from 
fits to the data for T=250 K, and calculated values from the 
shell model, which are the dashed lines. 



Energy spectra as a function of momentum and tem- 
perature were measured along the T-X direction. The 
phonon dispersion curves, obtained from fits to our data, 
are shown in fig.[3l^a). Some typical spectra can be found 
in the inset. We also overlay the dispersion curve calcu- 
lated using a shell model allowing for displacement of the 



atom's core and outer electron shell [32], implemented 
with the OpenPhonon code [33]. As shown in Fig. [3ja), 
the theoretical intensities of each mode, as represented 
by the thicknesses of curves, can vary greatly within the 
zone, to the extent that two modes may be mistaken for 
one when the intensity of one falls and the other rises. For 
example, modes in the 15-20 meV range which broaden at 
^=0.26 are shown in the inset of Fig.[3fa), and the broad- 
ening can be qualitatively explained by an anticrossing of 
two modes, predicted in the model. As further verifica- 
tion of the model, we plot in fig. [3fb) the ^-dependence 
of the relative integrated mode intensities at T=250 K, 
of the transverse acoustic mode and high-energy mode(s) 
that can be found in the 14-21 meV range, and normal- 
ize by the sum of the two intensities. By comparing the 
relative intensities of the modes, systematic error in the 
intensity that may arise from the beam hitting different 
parts of the sample as q is changed is canceled out. Also 
plotted is the corresponding calculated intensity fraction 
from the shell model for modes within the prescribed 
range, which is found to be in excellent agreement. We 
should note that our model parameters were originally 
optimized using SrTiOs dispersion data along the T-X 
and F-M directions [3J, |35|- The only subsequent mod- 
ification of the parameters was to change the mass (and 
atomic scattering factors affecting intensity) from Sr to 
Eu , whose most noticeable effect was a significant low- 
ering of the acoustic mode frequencies. 

The calculated eigenvectors describing the atomic 
motion are shown in table 1. The first two entries are for 
the lowest two optical modes calculated for g=0, TOi 
and TO2. The last entry in the table is the lowest energy 
mode at the R point. It consists of oxygen motions in 
the plane of the faces of the conventional cubic cell, and 
is a linear combination of degenerate modes of the form 
Oij=—Oji where Oij is the displacement of the oxygen 
on the z-facing plane in the direction. Such modes 
condense into the 14/ mem phase and correspond to the 
i?-point soft mode [19]. Table 1 confirms that the TOi 
mode is of "Slater" chatacter, similar to the soft mode in 
SrTiOa [35,^], while TO2 has Eu and Ti motions 180° 
out of phase. Looking at the temperature dependence 
of the modes in fig. O the energy of the TOi mode 
is generally the most sensitive to temperature in the zone. 

To investigate the temperature dependence of the 
TOi mode near the zone center, we changed to the (4 
1-1) and (5-1 -2) zones where our calculations predict 
a stronger intensity of that mode at the F-point. We 
note that even in these favorable zones, the mode is 
relatively weak compared to background coming from 
the tail of the acoustic modes, at high temperatures. 
Background-subtracted spectra, normalized for Bose 
factor and further by acoustic mode peak, are plotted in 
fig. [H for constant ^^=0.05 at selected temperatures. This 
latter normalization was to use the acoustic intensity 
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TABLE I: Table of calculated cubic-phase eigenvectors for 
the modes TOi at around 11.7 meV and TO2 at 20 meV at 
the F-point, and a q=(0.5 0.5 0.5) i?-point mode, denoted 
TAit:- The columns are the atom, the equilibrium position, 
followed by the displacements for each mode in units of lattice 
constant. 



Atom Position 



TOi 



TO2 



TAj 



Eu 
Ti 
O 
O 
O 



V 2 ' 2 ' 2 / 
(0,0,0) 

(io,o) 

(0,1,0) 
(0,0,1) 



(-0.007, 0, 0) 
(-0.011,0,0) 
(0.018,0,0) 
(0.038,0,0) 
(0.038,0,0) 



(-0.011, 0, 0) 
(0.03,0,0) 
(0.016,0,0) 
(0,0,0) 
(0,0,0) 



(0, 0, 0) 

(0,0,0) 
(0,-0.036,0.018) 
(0.036,0,-0.018) 
(-0.018,0.018,0) 



as a relatively stable reference against different sample 
positions at different temperatures, which could affect 
intensity. The ferroelectric mode energy vs. tempera- 
ture, measured over two zones, is plotted in the inset, 
including energies from additional spectra not plotted in 
the main figure. There is a general lowering of energy 
as the temperature is lowered, which is same order of 
magnitude previously observed in infrared reflectivity 
[rzirfisl (we estimate a '^4.5 meV shift by extrapolation 
to g=0, based on the dispersion in Fig. [3j), and the 
trend is approximately linear down to about 70-100 K 
after which the softening stops. The low temperature 
dispersion near ^=0 is plotted in fig. [3] as diamonds. 
Although high elastic background near zone center 
prevents measurement at exactly ^=0, the dispersion 
down to at least ^=0.02 is somewhat more linear than 
the room-temperature model would suggest for that 
mode in Fig. [3fa), where the dispersion instead flattens 
out as q is lowered. This linear dispersion characteristic 
of a soft mode [1|]. 

We observe no resolvable energy change below the 
magnetic transition TAr=5.7 K. This result may initially 
be surprising (not to mention disappointing), given the 
dramatic drop in dielectric constant measured in [sl, 
however we note that the magnitude of the dielectric 
constant drop within 1 K of the transition is less than 
0.3 %. A recent study of the optic mode using infrared 
reflectivity on thin EuTi03 films [37] films did in fact 
show an increase in phonon frequency below T/v, but in 
the vicinity of the transition, the change is less than 1%, 
which is indeed less than the present experiment could 
resolve for this mode. 

In summary, we prepared a single domain crystal of 
EuTiOa and observed it to undergo a transition into a 
structurally long-range ordered phase at r=287±l K, 
consistent with the 14/ mem structure from our powder 
refinements. The combined phonon and super-lattice 
Bragg peak intensity observations, which have consistent 
To, clearly prove a displacive structural transition. 
In this respect, the picture is considerably simpler 
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FIG. 4: (Color Online) Selected phonon spectra near the 
gamma point at g=0.05. Energy offsets were accounted for 
through fits of elastic intensity and acoustic modes on either 
side of huo^O where applicable. A smooth background was 
subtracted from the curves, whose intensities are corrected 
for Bose factor. The intensities were further normalized by 
the Bose-factor corrected transverse acoustic mode intensity. 
The Q position used was (5-1 -2.05) except for T=12 K which 
was taken at Q=(4 1 -1.05). Fitted energy vs. temperature is 
plotted in the inset for the (4 1 -1.05) and (5-1 -2.05) zones 
(circles and squares, respectively). 

than other recent diffraction studies in the literature 



jU, 23, 38|, where in some cases time-dependent incom- 
mensurabilities were observed. Sensitivity to sample 
preparation conditions may be a factor. We note that 
our data was taken over multiple beamtimes, and the 
phonons were basically consistent throughout. Near the 
structural transition, the i?-point soft phonon exhibits 
almost identical critical behavior to SrTiOa, with the 
same phonon energy vs. reduced temperature curve, 
suggesting a universal scaling relation for the soft mode 
energy. The lowest energy optical mode exhibits soft 
mode behavior consistent with displacive-type dynamics. 
A mass-corrected SrTiOa-based shell model can fairly 
accurately describe the frequencies and intensities of 
EuTiOa phonons, which further emphasizes the close 
similarity in the two systems. 
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